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SUMMARY 


A  series  of  simple  theoretical  models  are  discussed,  representing 
the  essential  features  of  flexible  aircraft  both  of  the  classical  and 
integrated  varieties,  in  trimmed  level  flight.  Assuming  linear  aero¬ 
dynamics  it  is  shown  that  the  maximum  speed  for  trimmed  flight  for  the 
classical  aircraft  is  determined  primarily  by  the  tail-plane  flexibility. 
For  the  integrated  aircraft  this  maximum  speed  occurs  when  both  the 
overall  aeroelastic  distortion  and  the  control  forces  become  very  large, 
the  effectiveness  of  the  control  to  trim  is  related  to  the  similarity 
of  the  weight  and  aerodynamic  lifting  distributions.  Detailed 
calculations  on  trimmed  plate  wings  of  varying  plan  form  further 
illustrate  these  general  points. 

The  aeroelastic  effects  of  these  trimmed  plate  wings  in  level  flight 
are  also  investigated  assuming  non-linear  aerodynamics.  It  is  shown 
that  in  general  there  are  two  positions  of  equilibrium  of  the  aircraft 
at  each  speed  although  it  is  possible  that  at  high  speeds  both  these 
positions  of  equilibrium  are  Imaginary.  The  stability  of  these  positions 
of  equilibrium  are  not  discussed  in  this  paper.  ^ 

SOMMAIRE  lXjk 

Ce  rapport  traite  d’une  slrie  de  noddles  thdoriwA^mples  represen¬ 
ted  les  caracteristiques  fondamentales  des  avions  souplV,  tant  classiques 
qu’  intrfgrds,  dans  le  vol  Iquilibrl  en  palier.  En  supposaimj*  alrodynamique 
Hn^alre,  11  est  ddmontrl  que  la  Vitesse  maximum  du  vol  4quiHbre  de 
1* avion  classique  est  dlterminle  principalement  par  la  souplease  de 
1* empennage.  En  ce  qui  concerne  1* appareil  int^gr^,  cette  vitesse 
maximum  se  produit  lorsque  la  distorsion  alrodlastique  globale  ainsi  que 
les  forces  des  gouvernes  deviennent  trds  Import antes:  1’ efficacltd  des 
gouvernes  en  vue  de  1*  equllibrage  est  fonctlon  de  la  similitude  des 
repartitions  du  polds  et  de  la  adrodynamique.  Des  calculs  detaillds 
relatifs  i  des  ailes  equillbrees  de  formes  en  plan  diverses  mettent 
plus  en  evidence  ces  points  gene'raux. 

Les  effets  aeroelastlques  de  ces  ailes  equilibrees  pour  le  vol  en 
palier  sont  egalement  etudlds  en  supposant  1'  adrodynaaique  nonlineaire. 

II  est  ddmontre  que  de  fa$on  generals  il  existe  deux  positions  d’equilibre 
pour  1*  avion  i  cheque  vitesse,  blen  qu’  il  soit  possible  qu'  £  des  vitesses 
elevens  ces  deux  positions  d’equilibre  soient  imaginaires.  La  stablllte 
de  ces  positions  d’equilibre  n’sst  pas  dlscutee  dsns  cette  communication. 


533.6.013.42 

JcSclal 


11 


CONTENTS 


Page 

SUMMARY  n 

LIST  OF  FIGURES  lv 

NOTATION  v 

1.  INTRODUCTION  1 

2.  MODEL  1.  REPRESENTING  THE  CLASSICAL 

AIRCRAFT  RITH  A  FLEXIBLE  MING  2 

3.  MODEL  2.  REPRESENTING  THE  CLASSICAL 

AIRCRAFT  RITH  A  FLEXIBLE  FUSELAGE  2 

4.  MODEL  3.  REPRESENTING  THE  CLASSICAL 

AIRCRAFT  RITH  A  FLEXIBLE  TAIL  3 

5.  MODEL  4.  REPRESENTING  AN  INTEGRATED  AIRCRAFT  4 

6.  THE  SLENDER  CONFIGURATIONS  6 

7.  SLENDER  DELTA  PLATE  RING  8 

8.  SLENDER  GOTHIC  PLATE  RING  11 

9.  CONCLUSIONS  11 

FIGURES  13 

DISCUSSION  A-i 

ADDENDUM:  Complete  List  of  Papers  la  Series 

DISTRIBUTION 


ill 


LIST  OP  FIGURES 


Page 

Pig. 1  Classical  aircraft  -  flexible  wing  13 
Pig. 2  Classical  aircraft  -  flexible  fuselage  13 
Pig. 3  Classical  aircraft  -  flexible  tail  14 
Pig. 4  Siaple  integrated  aircraft  14 
Pig. 5  Notation  for  slender  aircraft  15 
Pig. 6  Equilibrium  positions  of  delta  plate  -  non-linear  aerodynamics  15 
Pig. 7  Gothic  plate  wing  16 
Pig. 8  Equilibrium  shapes  of  gothic  plate  -  linear  aerodynamics  16 
Pig.  9  Equilibrium  shapes  of  gothic  plate  -  non-linear  aerodynamics  17 


ir 


SOME  STATIC  AEROELASTIC  CONSIDERATIONS  OP 
SLENDER  AIRCRAFT 

8.J.  Hancock* 


1.  INTRODUCTION 

The  techniques  for  the  calculation  of  aeroelastic  effects  have  been  formulated  and 
perfected  during  the  evolution  of  the  conventional,  or  classical  aircraft,  comprised 
of  distinct  wings,  fuselage  and  tail  and  acted  on  by  aerodynamic  forces  which  depend 
linearly  on  incidence.  We  are  now  in  a  state  of  transition  from  this  conventional 
subsonic  aircraft  to  the  supersonic  aircraft  of  the  future.  Hie  latter  is  expected 
to  differ  from  the  former  in  two  Important  respects.  First,  the  supersonic  con¬ 
figuration  will  be  integrated  into  a  single  unit  combining  the  roles  of  wings, 
fuselage  and  tail.  Second,  the  current  supersonic  aerodynamic  design  philosophy  of 
incorporating  flow  separations  and  the  formation  of  leading-edge  vortices  as  essential 
parts  of  the  flow  field  implies  that  the  aerodynamic  forces  will  be  non  linear 
functions  of  incidence.  Therefore,  from  an  aeroelastic  point  of  view  two  important 
questions  may  be  posed.  First,  which'  of  our  aeroelastic  concepts  based  on  background 
of  knowledge  built  up  with  the  experience  on  conventional  aircraft,  may  be  safely 
extrapolated  to  the  integrated  aircraft,  and  which  of  our  concepts  have  to  be  thought 
about  again  and  re-defined?  Secondly,  what  are  the  additional  aeroelastic  effects 
due  to  the  presence  of  the  non  linear  aerodynamic  forces?  In  this  Report  we  only 
begin  to  answer  these  questions  by  consideration  of  some  simple  static  aeroelastic 
effects. 

As  an  indication  of  the  reality  and  depth  of  these  questions,  a  difficulty  arises 
straight  away  over  definitions.  In  the  past,  static  aeroelasticity  has  been  defined 
in  terms  of  structural  and  aerodynamic  forces  only.  This  definition  tends  to  lose 
Its  meaning  when  applied  to  the  integrated  configuration,  so  the  following  definitions 
are  suggested.  Static  aeroelasticity  is  the  domain  of  aeroelasticity  in  which  all 
the  forces  (inertial,  structural  and  aerodynamic)  are  present  but  are  independent  of 
time.  Dynamic  aeroelasticity  is  the  domain  of  aeroelasticity  in  which  all  the  forces 
(inertial,  structural  and  aerodynamic)  are  present  but  are  dependent  on  time.  So 
static  aeroelasticity  refers  to  the  equilibrium  problem  of  a  steady  manoeuvre,  dynamic 
aeroelasticity  refers  to  stability  and  response  problems.  Static  stability  is  a 
Halting  case  of  dynamic  stability  and  is  not  included  under  the  present  definition 
of  static  aeroelasticity. 

In  this  Report  we  consider  only  static  aeroelasticity,  that  is,  the  equillbrius 
problem,  involving  a  steady  controlled  manoeuvre.  Since  one  manoeuvre  is  as  re¬ 
presentative  as  any  other  we  shall  lialt  our  discussion  to  the  triaaed  level  flight 
condition.  Since  we  are  Interested  in  relating  the  aeroelastic  effects  on  conventional 
aircraft  to  those  on  the  slender  Integrated  aircraft  we  shall  first  review  the  basic 
aeroelaatic  effects  on  conventional  aircraft  by  reference  to  some  extremely  simple 
models  which  idealise  the  various  flexibilities  of  this  type  of  aircraft.  Then  we 
shall  look  at  some  overall  effects  of  the  integrated  configurations,  again  hy  reference 
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to  simple  models  and  in  this  case  we  shall  include  a  consideration  of  the  non  linear 
aerodynamics. 


2.  MODEL  1.  REPRESENTING  THE  CLASSICAL 
AIRCRAFT  1ITH  A  FLEXIBLE  WING 

An  attempt  is  made  in  this  model  to  illustrate  some  of  the  main  ideas  regarding 
wing  distortion  on  a  conventional  aircraft  consisting  of  separate  wings,  fuselage  and 
tail  unit.  It  is  assumed  that  the  wings,  fuselage  and  tail  are  all  rigid,  but  that 
the  mounting  of  the  main  wings  to  the  fuselage  is  elastic,  allowing  the  wings  to  have 
the  single  degree  of  freedom  of  twist  only.  The  torsional  stiffness  of  the  wing  (of 
area  S„  )  relative  to  the  fuselage  is  denoted  by  (lb  ft/rad)  .  The  Incidence  of 
the  fuselage  is  denoted  by  a  .  The  angular  distortion  of  the  wing  relative  to  the 
fuselage  is  denoted  by  Q  ,  assuming  that  8=0  corresponds  to  the  unloaded 
equilibrium  attitude.  This  configuration  is  shown  in  Figure  1. 

For  the  purposes  of  this  model  it  is  assumed  that  C*0  is  zero.  This  assumption 
will  not  invalidate  the  qualitative  conclusions  derived  from  this  model.  Therefore  I 

in  steady  trismed  level  flight 


Wing  Lift  =  Lj 


Tail  Lift  =  4 


JipV2 S^a0  (a  +  8) 


(1) 


(2) 


Torsional  Distortion 


Ljec 

■s 


(3) 


Fuselage  Incidence  =  a 


)4/0V*  Sw 


_  !  HpV2  B||ec  a0 

ao  L  u0 


(4) 


(The  notation  is  given  in  Fig. 1).  From  the  expression  for  the  tall  lift  Lf  the 
elevator  angle  to  trim  can  be  calculated.  Thus  in  this  model  the  angular  distortion 
6  is  Independent  of  speed  but  the  fuselage  incidence  a  decreases  with  speed. 


Thus  from  a  trim  point  of  view  wing  distortion  is  not  a  critical  parameter. 

This  is  true  even  when  more  sophisticated  models,  representing  wing  flexibilities  in 
a  more  practical  manner,  are  assumed. 


3.  MODEL  2.  REPRESENTING  THE  CLASSICAL 
AIRCRAFT  HITS  A  FLEXIBLE  FUSELAGE 

In  this  model  we  look  at  the  flexibility  of  the  fuselage.  Again  to  illustrate  the 
main  points  we  consider  an  idealised  veraioo  of  the  problem.  The  fuselage  is 
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represented  by  an  elastic  beam  of  uniform  cross  section  and  uniform  weight  distribution 
(cr  lb/ft)  .  Hie  lifts  on  the  rigid  wing  and  tail  (L,  and  LT)  together  with  the  wing 
and  tail  weights  (Ww  and  ft)  act  at  the  ends  of  the  fuselage.  We  are  neglecting  the 
part  of  the  fuselage  ahead  of  the  wings  and  aft  of  the  tail.  The  wing  section  is 
assumed  to  be  symmetrical  (l.e.  angle  of  zero  lift  is  zero)  with  its  centre  line 
parallel  to  the  local  fuselage  centre  line,  i.e.  the  wing  incidence  is  the  same  as 
the  incidence  of  the  front  of  the  fuselage.  This  is  shown  in  Figure  2. 

The  steady  trimmed  level  flight  condition  gives 
L,  =  »,+%  crl  =  KpV2  Swa0aXB0 


h  =  WT+fc<^  =  [a,  (axal  -  ^ax»o)  +  a^] 

The  fuselage  distortion  f(x)  may  be  found  by  solving  the  differential  equation 
for  the  loaded  beam  with  the  appropriate  end  conditions.  The  main  results  however  may 
be  deduced  directly  from  physical  conditions.  These  are  (i)  since  L*  is  constant 
then  (a)  ,0  will  be  the  same  whether  the  beam  fuselage  is  rigid  or  flexible,  it  is 
independent  of  the  beam  flexibility;  (ii)  since  the  loading  in  the  trimmed  state  is 
independent  of  speed  the  fuselage  distortion,  measured  relative  to  the  axis  (a)z„0  , 
is  Independent  of  speed,  lhus  from  a  trim  point  of  view  fuselage  flexibility  is  not 
a  critical  parameter. 


(5) 

(6) 


4.  MODEL  3.  REPRESENTING  TRE  CLASSICAL 
AIRCRAFT  MITH  A  FLEXIBLE  TAIL 

10  coa«>lete  our  analysis  of  the  conventional  configuration  we  should  discuss  the 
distortion  characteristics  of  a  flexible  tail  plane.  For  this  model  we  take  a  rigid 
wing  (symmetrical  section,  zero  angle  relative  to  fuselage),  rigid  fuselage  and  a 
rigid  tail  plane  with  a  flexible  mounting  of  torsional  stiffness  U6.  (lb  ft/rad) 
between  the  tail  plane  and  fuselage.  This  is  shown  in  Figure  3. 


In  steady  trimmed  level  flight 


^w  =  * 


1,  +  h 


=  )i/0V2Swa0a 


(7) 


=  1 


Iw  +  lT 


=  #/0V*ST  (ml  [a  +  d  -  e]  +8,77) 


The  structural  equilibrium  of  the  tall  is  given  by 

Ltj, ec  -  =  me6 


(8) 


(8) 
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where  is  the  moment  about  the  tail  aerodynamic  centre  (+  nose  down)  due  to  the 
elevator  deflection  17  defined  by 

Atp  =  J4/oV2STcTmT7  (10) 


The  effect  of  MT  on  the  overall  equilibrium  conditions  is  usually  neglected.  The 
substitution  of  equations  (9)  and  (10)  into  Equation  (8)  gives  the  elevator  angle  to 
trim. 

[i^  -  )4/oV2STa1a  (1  -  de/3a)]  [m5  -  \)4pV2 ST cTaj  m/a2] 

71  )4pV2  SjBj  [ug  -  )4pV2  STcTa1m/a2] 

where 

K  =  ea*  h _ 

m  1^  -  )4pV2STa1a(l  -  3e/3a) 


The  elevator  angle  to  trim  tends  to  infinity  as  V  -  Vc  where 


c  ^)4pSTcTa1  m/a2 

The  speed  Vc  is  therefore  the  maximum  speed  theoretically  possible  for  trimmed  level 
flight.  Alternatively,  if  the  elevator  effectiveness  for  trim  is  defined  as  the 
elevator  required  to  trim  (l.e.  Equation  11)  divided  by  the  elevator  required  to  trim 
with  a  completely  rigid  tail  (i.e.  Equation  li  with  m0  infinite),  then  Vc  is  the 
speed  at  which  there  is  no  control  effectiveness  to  trim.  This  is  sometimes  referred 
to  as  the  ‘elevator  reversal  speed’  but  this  definition  is  not  really  admissible  since 
the  aircraft  cannot  attain  the  speed  at  which  it  cannot  be  controlled,  and  so  the 
concept  of  reversal  of  control,  which  implies  flight  through  the  critical  region, 
is  not  relevant. 

Note  that  when  A.  =  1  there  is  no  loss  of  control  effectiveness  except  at  the 
critical  speed  VQ  .  Therefore  from  a  trim  point  of  view  .tail  flexibility  is  the 
critical  parameter  which  determines  the  overall  aeroelastlc  limiting  conditions  for 
this  controlled  condition  of  flight. 


(12) 


S.  MODEL  4.  IEPME8ENTIN8  AN  INTE6EATED  AIBCBAFT 

Ve  now  consider  an  idealized  integrated  system  in  which  the  wing,  fuselage  and 
tail  are  combined  together,  as  they  weld  be  for  a  slender  configuration  with  the 
aerodynamic  loading  distributed  over  the  whole  of  the  length.  Ve  take  as  our  model 
a  staple  beam  of  uniform  cross-section  (uniform  bending  stiffness  distribution  El, 
uniform  weight  distribution  a  lb/ft,  total  length  l  )  under  a  simplified  linear 
aerodynamic  load  distribution  which  is  assumed  to  be  KpV2a(x)  ,  where  a(x)  is 
the  local  incidence  and  K  is  taken  to  be  constant.  To  simulate  the  force  on  the 
elevator  control,  a  force  P  is  applied  at  the  trailing  edge  of  the  system,  ensuring 


that  the  overall  equilibrium  conditions  (Lift  =  Weight,  Moment  =  0)  are  satisfied. 
This  system  is  shown  in  Figure  4. 
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In  steady  trimmed  flight  the  distortion  function  y(x)  satisfies  the  differential 
equation 


El  ylT  =  -or  -  KpV* y' 

The  boundary  conditions  are 

y"'lx-0  =  y"lx-o  =  ylx«i  =  0 


(13) 


(14) 


'  xm  l  ~ 


P 

ii 


(15) 


It  can  be  shown  that  any  solution  of  Equation  (13)  which  satisfies  the  boundary 
conditions  (14)  and  (15)  automatically  satisfies  the  overall  equilibrium  conditions. 


The  solution  of  Equation  (13)  is 


-y'  =  a(x)  = 


K  pVz 


+  Ae*x  +  e"^x  ^B  cos  —  Ax  +  C  sin  Ax^ 


(16) 


where 

.  -S £* 

El 


(17) 


and  A  ,  B  and  C  are  arbitrary  constants.  The  substitution  of  Equation  (16)  into 
the  three  boundary  c (Hid it  ions  (14)  gives  three  simultaneous  equations  for  A  ,  B  and 
C  (the  coefficients  of  which  depend  on  A  ).  Thus  the  general  solution  is 

A  =  B  =  C  =  0 

except  at  a  critical  velocity  Vc  ,  corresponding  to  the  value  of  A  when  the 
determinant  of  the  simultaneous  equations  is  zero.  At  this  speed  A  ,  B  and  C  are 
indeterminate,  although  B  and  C  can  be  found  in  terms  of  A  . 

Equation  (15)  gives  the  control  force  P  once  the  distortion  is  calculated.  Thus 
in  general  P  is  zero  except  at  V  =  V  where  P  depends  on  the  coefficient  A 
(assuming  that  B  and  C  have  been  calculated  in  terns  of  A  ). 

Therefore,  for  this  model  in  the  trimmed  state  there  is  no  distortion  until  a 
critical  speed  is  reached  when  a  distortion  can  exist,  depending  on  the  size  of  an 
arbitrary  control  force  P  .  The  absence  of  distortion  for  speeds  below  Vc  is 
explained  physically  in  the  present  model  by  the  fact  that  the  lift  distribution 
exactly  balances  the  weight  distribution.  In  order  to  deduce  further  information 
from  the  model  it  is  neoessary  to  alter  this  balance  of  aerodynamic  and  weight 
loadings.  This  nay  be  done  by  modifying  the  uniform  weight  distribution  in  this  model 
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(as  shown  in  Pig. 4)  to  an  asymmetric  weight  distribution.  If  it  is  now  assumed  that 
the  weight  distribution  is  taken  to  be 


cr 


then  the  distortion  function  y(x)  satisfies  the  differential  equation 


(18) 


El  yiv 


-  KpV2  y ' 


(19) 


The  boundary  conditions  are  the  same  as  before,  l.e.  Equations  (14)  and  (15).  The 
solution  of  Equation  (19)  is 


-y' 


a(x)  = 


o-q  (2  -  x/l) 
Kp  V2 


+  AeXx  +  e 


-KXx 


(  A  v5  \ 

[  B  cos  —  Ax  +  C  sin  —  Ax ) 

V  2  2  J 


(20) 


where  X  is  the  same  as  before  (i.e.  Equation  17).  The  substitution  of  Equation  (20) 
into  the  boundary  conditions  (14)  again  gives  three  simultaneous  equations  for  A  , 

B  and  C  which  in  this  case  gives  finite  non-zero  values  for  these  coefficients  for 
each  value  of  speed  (or  X  ).  The  control  force  P  can  then  be  determined  by  sub¬ 
stituting  these  known  values  of  A  ,  B  and  C  into  Equation  (20)  and  then  satisfying 
condition  (15).  However,  A  ,  B  and  C  Increase  as  X  approaches  the  critical 
value  at  which  the  determinant  of  coefficients  of  A  ,  B  and  C  in  the  simultaneous 
equation  is  zero.  At  the  critical  speed  Vc  the  distortions  tend  to  infinity  as  the 
control  force  P  becomes  very  large. 


It  is  noted  that  the  critical  speed  is  the  same  for  both  the  weight  distributions. 
Therefore  the  critical  speed  is  independent  of  the  weight  distribution  but  the  control 
effectiveness  to  trim  below  this  critical  maximum  trim  speed  depends  on  the  out-of- 
balance  of  the  aerodynamic  lifting  and  weight  distributions. 


It  could  be  argued  that  the  aeroelastlc  effects  associated  with  the  tail  flexibility 
on  conventional  aircraft  are  similar  to  the  overall  effects  of  the  whole  of  the 
Integrated  slender  aircraft. 


8.  THE  SLENDER  CONFISCATIONS 

The  previous  simplified  models  have  given  a  qualitative  idea  of  the  main  aero- 
elastic  effects  which  we  night  expect  to  occur.  At  this  stage  therefore  it  is 
necessary  to  extend  these  models,  especially  those  representing  the  integrated 
configurations,  to  Incorporate  the  particular  features  of  the  slender  aircraft,  both 
structurally  and  aerodynamically.  The  slender  configuration  is  shown  in  Pigure  5, 
together  with  the  associated  notation. 

It  is  assumed  that  this  aircraft  will  bend  as  a  beam.  For  a  given  structure  the 
bending  stiffness  distribution  EI(x)  and  the  lengthwise  weight  distribution  w(x) 
are  assumed  to  be  known.  The  main  difficulty  at  this  stage  is  the  general  lack  of 
knowledge  of  the  aerodynamic  loading  on  streamwise  cambered  slender  wings.  However, 
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It  is  assumed  in  this  Report  thst  the  loading  distribution  l(x)  is  given  by  the 
equation 


l(x) 


7T8*(x)a(x)  +  4a(x)|a(x)| 


dx 


(21) 


for  the  general  shape  wing  (l.e.  s(x)  and  incidence  distribution  a(x)  ).  It  is 
assumed  here  that  a  =  0  is  the  zero  lift  case.  This  is  an  Intuitive  approximation 
based  on  the  following  observations: 


(1)  If  the  non-linear  terms  are  absent  (i.e.  if  there  are  no  separation  and 
associated  vortices)  then  the  loading  is  given  by  linearized  slender-wlng 
theory  which  states  that 

l(x)  =  pV2^-  [7rs*(x)a(x)] 
dx 

for  general  distributions  s(x)  and  a(x)  . 

(il)  The  only  reliable  formula  incorporating  the  leading-edge  vortices,  based  on 
both  theoretical  and  experimental  grounds,  is  for  plate  (or  symmetrical) 
delta  aerofoils.  This  formula  is 

l(x)  =  pV2  [27rk2  xa  +  4  kxa|a|  ] 

where  s(x)  =  kx  and  a(x)  =  a  which  is  constant. 

Equation  (21)  is  correct  for  these  two  particular  cases  and  it  is  hoped  that  its 
generalized  fora  will  give  a  use  Ail  loading  formula  which  represents  the  essential 
loading  characteristics  of  deformed  slender  shapes. 

The  Overall  equilibrium  conditions  which  must  be  satisfied  in  the  trimmed  state 
are 

jT  [  I(x)  -  w(x)]  dx  +  P  =  0  (22) 

and 

J  [l(x)  -  w(x)]  [x  -  c]  dx  =  0  (28) 

The  differential  equation  for  the  deformed  aircraft  shape  Mx)  is 

dx* 


(II  i")  =  *pV*— 

uX 


7TB 


+  4i 


•x 

I'li'lj  s(x)  dx 

Jo 


-  w(x) 


(84) 


with  the  boundary  conditions 
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biA"|x.0  =  BiA"|x. 


-  (BI*")|X.0  =  0 


(25) 


-  <EIA")|X.0  =  -P 


(26) 


This  Equation  (24)  nay  be  regarded  as  a  third-order  equation  for  the  incidence 
distribution  a(x)  satisfying  the  three  boundary  equations  (25).  From  this  solution 
the  control  force  P  for  overall  equilibriua  is  then  calculated  from  Equation  (26). 
On  integration,  applying  the  appropriate  boundary  conditions  Equation  (24)  becomes 


El  a'  (g)  =  -pV2  c*1 


1 J ^  ^*2(f)  otfj)  +  4a(^j)  |a(^)  |  d£ 


*x 


(27) 


3  ff  P1 

+  c*  ' 
Jo  Jo 


*(f) 4ft df 


where  £  =  x/c 

with  the  boundary  condition 


BIci'  |  =0  (28) 

Equation  (27)  is  completely  general,  and  it  can  be  solved  for  any  distribution 
■I  .  *(£)  .  s(£)  etc.  by  a  simple  step  by  step  procedure  fros  f  =  0  to  f  =  1  . 

However,  to  Illustrate  the  sain  qualitative  trends  fros  this  equation  we  shall 
concentrate  once  again  on  sons  staple  slender  nodels.  These  are  taken  to  be  slender 
flat  plates,  first  a  slvle  delta  and  the  second  a  gothic  outline.  Although  these 
plates  do  not  have  the  type  of  distributions  of  weight  and  stiffness  which  are 
expected  in  a  typical  supersonic  project  they  are,  at  least,  consistent  in  their 
relationship  between  stiffness  and  weight  sad  they  will  show  the  sain  seroelsstie 
trends  of  distortion  in  trlnaed  flight. 


7.  8LENDEK  DELTA  PLATE  WIN6 

We  now  consider  the  siaple  case  of  a  slender  delta  wing  of  constant  thickness  t 
and  specif ic  weight  cr  per  net  area  (lb/ft2)  trianed  for  level  flight  hy  a  control 
force  P  at  the  trailing  edge.  In  this  case  therefore 

s(x)  =  Kx  =  Eof  (28) 

The  substitution  of  the  appropriate  stiffness  and  weight  distributions  into  the 
differential  Equation  (28)  gives 
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a  .  -  -2>(£ifXr-3 


A,  = 


The  boundary  condition  (28)  is 


a'(l)  =  0 


We  may  consider  this  problem  in  two  parts,  namely  with  and  without  the  non-linear 
part  of  the  aerodynamic  loading,  that  is,  considering  the  problem  in  the  presence  of, 
and  in  the  absence  of,  the  leading  edge  vortices. 

(i)  Linear  Aerodynamics 

The  differential  equation  in  this  case  is 

«•’  »  <»> 


and  <%(£)  must  satisfy  the  end  condition  given  by  Equation  (31).  The  first  obvious 
solution  which  satisfies  the  boundary  condition  is  that  a  is  constant  and  is  given 
by 


This  solution,  however,  is  not  unique.  This  is  shown  by  obtaining  the  general  series 
solution  of  Equation  (32),  which  is 

a  -  a0  -  (A<x0-  At)^#f  -  jl<#t  +  3»  #*  9*J 

Boundary  condition  (81)  gives 

f  /i  A  A2  \ 

0  <M> 

Therefore  a  is  constant,  given  by  Equation  (33),  except  at  the  critical  speed 
(or  A  )  when  the  second  bracket  in  Equation  (34)  is  sero,  and  at  this  speed  a0  is 
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arbitrary.  This  implies  that  the  control  force  P  is  zero  except  at  the  critical 
speed  where  any  control  force  P  nay  be  applied,  keeping  the  plate  in  trin. 

This  problen  is  therefore  equivalent  to  Model  4,  since  physically  the  weight  and 
aerodynamic  lifting  distributions  again  exactly  balance.  The  critical  speed  from 
Equation  (34)  is  given  by 


A  =  33 


(35) 


(li)  Non-Linear  Aerodynamics 

The  differential  equation  in  this  case  is  the  tall  Equation  (30),  together  with 
the  boundary  condition  (31).  Again  in  this  case  a  trivial  solution  exists  giving  a 
constant,  where  a  is  given  by 


a 


±1 

A 


(36) 


Physically  this  is  because  the  presence  of  the  vortices,  according  to  the  present 
aerodynamic  assumption  (l.e.  Equation  21),  do  not  alter  the  distribution  of  loading; 
they  only  increase  the  intensity  of  loading.  So  the  plate  takes  up  a  smaller  angle 
of  incidence,  without  distortion.  And  again  solution  (36)  is  not  unique.  A  series 
expansion  can  be  applied,  assuming  a  >  0  for  0  4  £  £  1  ,  giving  a  solution  of  the 
following  form 

The  boundary  condition  (31)  gives  the  result  that 


P  (f  A.  Alt  a0)  (37) 


Therefore, at  each  speed,  in  addition  to  the  solution  (36),  there  is  the  possibility 
of  another  distortion  shape  if 


F'  (1,  A,  Alt  a0)  =  0  (39) 

There  seems  to  be  at  least  one  simple  solution  of  Equation  (39)  so  that  in  general  at 
each  speed  there  appear  to  be  two  distinct  positions  of  equilibrium.  As  the  speed 
increases  these  two  distortion  shapes  merge  together  until  a  certain  speed  is  reached 
at  which  both  distortion  shapes  coincide.  Above  this  speed  the  two  distortion  shapes 
diverge  from  each  other.  This  speed  at  which  the  two  shapes  coincide  might  be 
regarded  as  a  critical  speed,  near  which  the  aircraft  would  tend  to  Jump  from  one 
position  of  equilibrium  to  the  other,  assuming  that  both  positions  of  equilibrium  are 
stable.  It  could  be  argued  that,  in  the  region  of  this  critical  speed,  if  one 
position  of  equilibrium  is  stable  and  the  two  positions  gradually  coincide  with 
increase  of  speed  the  other  position  of  equilibrium  is  stable.  This  question  of 
stability  cannot  be  resolved  at  this  stage  and  requires  further  investigation  before 
a  definite  conclusion  may  be  readied  about  the  relationship  of  these  positions  of 
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equilibrium.  In  this  Report  we  have  only  established  their  existence.  The  critical 
velocity  of  this  plate,  froa  Equation  (38)  cones  out  to  be 

A  =  19.1  (40) 

It  is  interesting  to  note  the  low  value  of  this  speed  compared  to  the  critical  speed 
in  the  absence  of  the  leading-edge  vortices  (Equation  35). 

The  equilibrlun  shapes  at  A  =  16.5  and  A  =  22  are  shown  in  Pigure  6. 


8.  SLENDER  OOTHIC  PLATE  UNO 

This  gothic  plate  wing  is  shown  in  Pigure  7.  Again  this  problem  was  Investigated 
with  and  without  the  non-linear  aerodynamic  terns;  in  this  problen  the  step  by  step 
numerical  techniques,  nentioned  earlier,  had  to  be  applied.  The  results  are 
sunaarized  below. 

(i)  Linear  Aerodynamics  H 

The  naxlnun  trln  speed  in  this  case  is  given  by 

A  =  44.5  (41) 

The  variation  of  the  equilibrlun  shapen  of  the  wing  at  values  of  below  the  critical, 
together  with  the  variation  of  control  force,  are  shown  in  Pigure  8. 

(11)  Non-Linear  Aerodynamics 

It  vis  found  in  this  case  that  at  low  values  of  A  two  positions  of  equilibrlun 
are  apparently  possible.  This  in  a  similar  result  to  the  delta  plate.  At  a  critical 
speed,  given  fay 


A  =  21.6  (42) 

both  positions  of  equilibrlun  coincide  to  a  single  position.  Above  this  critical 
speed  no  positions  of  equilibrium  are  apparently  possible.  Shapes  of  equilibrium  are 
shown  in  Pigure  9. 


9.  CONCLUSIONS 

The  sain  conclusions  regarding  the  aeroelastic  effects  in  trlswed  level  flight  nay 
be  sussarised  as  follows; 

(1)  On  a  conventional  aircraft  it  is  prisarily  the  tall  flexibility  which  is 
responsible  for  the  loss  of  control  effectiveness  to  trls  resulting  in  a 
saxlsun  tris  speed. 

(ii)  lbs  assumption  of  linear  aerodynaalos  on  an  integrated  aircraft  results  in 
a  loss  of  control  effectiveness  to  tris  with  increasing  speed,  giving  a 
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maximum  trim  speed  (1. e.  zero  control  effectiveness)  at  which  the  dis¬ 
tortion  of  the  whole  airfraae  Is  very  large.  This  Uniting  or  critical 
speed  is  independent  of  the  weight  distribution,  but  the  control  effectiveness 
effectiveness  to  trim  depends  on  the  out -of -balance  of  the  aerodynamic  and 
weight  distributions. 

(ill)  Proa  the  assumption  of  non-linear  aerodynamics  on  an  integrated  aircraft  it 
appears  that  at  low  speeds  two  independent  positions  of  equilibrium  exist. 
These  two  modes  of  distortion  coincide  at  a  certain  critical  speed.  Above 
this  critical  speed  there  may  be  either  two  positions  of  equilibrium  or  no 
positions  of  equilibrium.  The  nature  of  this  critical  speed  has  not  been 
satisfactorily  explained  in  this  Report;  this  can  only  be  done  by  reference 
to  the  stability  of  the  various  modes  of  equilibrium.  This  investigation 
is  at  present  under  consideration. 

It  is  interesting  to  note  that  the  critical  speed  with  non-linear  aerodynamics  is 
considerably  less  than  the  critical  speed  with  linear  aerodynamics.  This  may  be 
relevant  to  the  current  ideas  on  the  performance  of  slender  aircraft  when  the  cruise 
condition  is  designed  to  have  little  or  no  flow  separation  at  the  leading  edges  so 
that  the  cruise  is  determined  by  linear  aerodynamic  considerations,  whereas  the  off- 
cruise  condition  will  have  to  contend  with  the  effect  of  leading -edge  vortices,  that 
is  the  effects  of  non-linear  aerodynamics. 

Finally,  this  paper  is  an  introductory  attempt  to  understand  some  of  the  fundamental 
principles  of  aeroelasticlty  relating  to  the  integrated  aircraft.  The  approach 
advocated  in  this  paper,  of  studying  simple  models,  produces  the  general  background 
which  is  desirable  before  any  detailed  design  work  is  entered  into. 
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Pig. 1  Classical  aircraft  -  flexible  wing 


Pig.  2  Classical  aircraft  -  flexible  fuselage 
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Pig.  3  Classical  aircraft  -  flexible  tall 
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PK.  4  Staple  integrated  aircraft 
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Pig. 5  Notation  for  slender  aircraft 
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equilibria  positions  of  delta  plate  -  non-linear  aerofrasales 


DISCUSSION 


G.H.  Lee  (U.K.):  Mr.  Hancock  has  shown  that  with  linear  aerodynamics  a  flat  plate  of 
gothic  planform  can  have,  in  general,  two  different  equilibrium  shapes,  though  at  a 
certain  speed  these  two  shapes  are  the  same.  However,  with  non-linear  aerodynamics 
no  solution  is  possible  at  speeds  beyond  that  for  which  the  two  shapes  coincide; 
he  therefore  calls  this  a  critical  speed. 

Since  the  difference  in  these  two  cases  is  due  to  the  incorporation  of  non-linear 
terms,  I  should  like  to  know  if  Mr.  Hancock  has  any  idea  as  to  how  sensitive  this 
speed  is  to  the  type  of  non-linear  terms  incorporated.  In  other  words,  would  you  get 
similar  critical  speeds  with  different  but  reasonable  non-linear  theories,  or  is  the 
critical  speed  going  to  be  very  sensitive  to  the  type  of  theory  used? 

Reply  by  Author:  Pirst,  I  must  correct  the  misapprehension  gained  by  Mr.  Lee  from 
my  lecture;  ■  with  linear  aerodynamics  there  is  only  one  position  of  equilibrium  at 
each  speed,  it  is  only  with  non-linear  aerodynamics  that  there  arises  the  possibility 
of  more  than  one  position  of  equilibrium.  However,  ay  paper  is  Intended  to  be  rally  a 
qualitative  introduction  to  the  type  of  aeroelastic  phenomena  which  nay  arise  on 
slender  aircraft  In  the  presence  of  non-linear  aerodynamic  forces;  such  critical 
speeds  will  occur  with  other  assumptions  for  the  non-linear  aerodynamic  loading 
although  it  is  premature  at  this  stage  to  give  reliable  quantitative  assessments  of 
these  critical  speeds.  It  should  be  remembered  that  quantitative  results  depend  not 
only  on  a  more  reliable  knowledge  of  the  non-linear  aerodynamics  but  also  on  more 
realistic  estimates  of  the  structural  and  inertial  loadings. 


H.H.B.M.  Thomas  (U.K.):  Regarding  Mr.  Lee’s  question  to  Dr.  Hancock,  the  point  to 
note  is  that  linear  aerodynamics  leads  to  the  cruise  as  the  critical  condition 
aeroelastically  with  the  flight  path  plan  proposed  for  such  aircraft.  The  query 
raised  by  the  present  paper  is  whether  off -design  condition  at  larger  Incidence  nay 
not  be  more  critical  because  of  non-linear  aerodynamics.  I  an  sure  Dr.  Hancock  will 
agree  that  his  papsr  points  to  the  existence  of  a  possible  problem  rather  than  an 
assessment  of  the  severity  of  the  problem  in  an  actual  design. 
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